Non-perturbative singularity structure of gluon and quark propagators 
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A gauge invariant, non-local observable is constructed in lattice pure gauge theory, which is 
identical to the gluon propagator in a particular gauge. The transfer matrix formalism is used to 
show that this correlator decays exponentially with eigenvalues of the Hamiltonian. This implies 
a gauge invariant singularity structure of the propagator in momentum space, permitting a non- 
perturbative definition of a parton mass. The relation to gauge fixing and the extension to matter 
fields are discussed. 



PACS: 11.15, 11.15.Ha, 12.38.Aw, 12.38. Gc 

The confinement problem of QCD consists of the dy- 
namical relation between perturbative parton physics at 
short distances and non-perturbative hadron physics at 
large distances. To understand the transition of a system 
from one regime to the other requires a non-perturbative 
study of the dynamics of colour degrees of freedom, which 
is encoded in the Green functions of quark and gluon 
fields, the fundamental degrees of freedom of the ac- 
tion. This is particularly essential in the context of finite 
temperature and density physics probed in experimen- 
tal heavy ion collisions, where creation of a "deconfincd" 
quark gluon plasma is expected, whose collective physical 
properties should be determined by parton dynamics. 

In perturbation theory one fixes a gauge and studies 
parton interactions directly. Because of their gauge de- 
pendence, field propagators are not physical observables. 
Nevertheless, physical information about the parton dy- 
namics is carried by their singularity structure. For ex- 
ample, the pole mass defined from the quark propagator 
is gauge independent and infrared finite to every finite 
order in perturbation theory [jlj. A similar result holds 
for the gluon propagator, provided an appropriate resum- 
mation of infrared sensitive diagrams has been performed 
. Gauge invariant resummation schemes have been de- 
signed to self-consistently compute the pole of the gluon 
propagator in three dimensions which is related to 
the "magnetic mass" regulating the non-abelian thermal 
infrared problem In a Hamiltonian analysis of the 
three dimensional gauge theory a gauge invariant com- 
posite gluon variable has been constructed, which in the 
weak and strong coupling limits yields a gluon mass gap 
as the lowest eigenvalue of the Hamiltonian (s| . 

However, perturbation theory and resummation meth- 
ods are limited by the requirement of weak coupling, and 
nothing is known about the non-perturbative existence of 
field propagator poles in a confining regime. Hence a non- 
perturbative analysis in the framework of lattice gauge 
theory is warranted. Unfortunately, numerical gauge fix- 
ing on the lattice || is plagued by several problems: 
First, it is difficult to fix a gauge uniquely and avoid 
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the problem of Gribov copies Q . Second, most complete 
gauge fixings (e.g. the Landau gauge) violate the posi- 
tivity of the transfer matrix, thus obstructing a quantum 
mechanical interpretation of the results. Because of these 
problems, many results from gauge fixed simulations have 
remained controversial. An overview with references to 
recent numerical work may be found in Q . 

For the extraction of quark masses gauge fixing can be 
circumvented by the methods of non-perturbative renor- 
malization, which make use of manifestly gauge invariant 
quantities |j) . In this approach quark masses are defined 
through PCAC relations rather than by a field correlator, 
but this does not permit to address the same question for 
gluons. 

In this letter, the transfer matrix formalism Jl(]Jll| is 
used to show that field correlation functions carry gauge 
invariant information non-perturbativcly. To this end, 
SU{N) pure gauge theory with Wilson action is consid- 
ered in d + 1 dimensions on a (aL) d x aN t lattice with 
periodic boundary conditions. 

A gauge invariant composite and local field carrying 
the quantum numbers of a gluon can be constructed in 
theories with a complex scalar iV-plet in the fundamen- 
tal representation, such as the electroweak sector of the 
standard model. In the pure gauge theory no such field 
is present. However, it is possible to construct complex 
functions of the gauge field f a [U],(a = 1,...N), with 
the same transformation behaviour from the covariant 
Laplacian, 
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The latter is a hermitian operator with a strictly positive 
spectrum, whose eigenvectors have the desired transfor- 
mation property f^ B (x) = g{x)f^ n \x). They provide a 
unique mapping U — > f[U] except when eigenvalues are 
degenerate or |/| = 0. In practical simulations the proba- 
bility of generating such configurations is essentially zero 
|p^| . These properties have been used previously for 
gauge fixing without Gribov copies [fl3[ and to construct 
blockspins for the derivation of effective theories MM . 



The eigenvectors are non-local in the sense that they 
depend on all link variables. In order to maintain the 
transfer matrix formalism the f(x) have to be local in 
time. This is achieved by considering the spatial Lapla- 
cian Dj{Ui] in Eq. ([!]), which then is defined in every 
timeslice individually and independent of Uq. 

The eigenvectors are used to construct an TV x N ma- 
trix Q(x) G SU(N) following Since Eq. (@) only 
determines them up to a phase, this leaves a remain- 
ing freedom in fl(x). In the case of SU(2), all eigen- 
values are two-fold degenerate due to charge conjuga- 
tion, and the two vectors to the lowest eigenvalue are 
combined into f2, which then is determined up to a 
global SU(2) rotation h. For SU(3) there is no de- 
generacy of the eigenvalues in general. In this case 
one solves for the three lowest eigenvectors to construct 
the matrix Q, which is then determined up to a factor 

h = diag(exp(iwi),exp(iw2),exp(iu;3)),X^^ = 0. This 
may be summarized by the transformation behaviour 



n 9 { x ) = ff (x)o(x)/i t (t), 



(2) 



where h(t) is free and may be different in every timeslice. 
We can now define composite link and gluon fields 



V ll (x)=a*{x)U ll (x)il(x + (i), 



(3) 
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(V^x) - V£{x)) - -Tr (V,(x) - Vfcx)) 



both transforming as Of(x) = h(t)Oi(x)hf (t), whereas 
V 9 (x) = h(t)V (x)h\t + 1). Hence the A, are gauge 
invariant under spatial transformations <?(x), but trans- 
form under time-dependent rotations corresponding to 
the residual symmetry of the spatial Laplacian. 

With Vo(x;ti,t2) denoting the temporal Wilson line 
connecting the sites (x, t\), (x, £2)1 we can now construct 
the manifestly gauge invariant operator 

0[U] = Tr [A 4 (x,0)K,(y;0,t)A i (x,t)F t (z;0,t)] . (4) 

It represents a correlator in t of the composite field Aj, 
where Wilson lines are inserted to ensure full gauge in- 
variance. Note that these may be placed at any y, z. 

In order to obtain a spectral decomposition of {0[U}) 
it has to be converted into a trace over quantum me- 
chanical states. In a Hilbert space formulation ]l5| , ^0| , pl| 
a spatial sublattice L d at a fixed time is considered, with 
link variables U (x, i). The wave functions form a Hilbert 
space TLq of all complex, square integrable functions 
ip[U] defined on the gauge group G: H = [L 2 (G)] dL " . 
Wave functions of physical states are gauge invariant, 
ip[U 9 } = i/j[U], forming a subspace H cHo- Any wave 
function ij>[U] £ Tlo can be projected on the physical 
subspace by means of the projection operator P , 



(5) 



The dynamics is introduced by means of the transfer ma- 
trix T, which translates wave functions by one lattice 
spacing in time. It is an integral operator 

iipt+x] = (iV) [U t ] 

= J Y[ d Ut(x,i) K[U t+u UMU t ] , (6) 



with kernel 



K[U, U']= JYl dW(x) exp -S t [U, W, U 1 ] 



(7) 



Here St is the action of two neighbouring timeslices, and 
the field W(x) is to be identified with the timelike links 
. Using the gauge invariance of St , one shows that 



T = PTn 



(8) 



where To is an integral operator with kernel Kq [U, U'] — 
exp—St[U,l,U'], which is invariant under time inde- 
pendent gauge transformations ,g(x). T is a bounded, 
sclf-adjoint operator with a strictly positive spectrum 
ifTTf , thus allowing to define a Hamiltonian jll],[l0| Hq = 
— 1 1 a In To, acting on the Hilbert space 7io- Through 
the projection T is defined on the gauge invariant sub- 
space 7i, on which there is a corresponding Hamiltonian 
H. Generally, the spectrum of H consists of the phys- 
ical particle states of the theory, which couple to local 
gauge invariant operators. In addition to these states, 
Ho contains also the spectrum of gauge field excitations 
in the presence of static sources, such as the static po- 
tential. The corresponding pure gauge wave functions in 
one timeslice have non-trivial transformation behaviour. 

We now define the multiplication operators 
{&affll>)[U] = n aP [U]il>[U] and (V(x,^)[[/] = 
F(x, i)tp[U], so that A(x,i) is a multiplication opera- 
tor as well. 

Since 0[U] is manifestly gauge invariant, its expecta- 
tion value is the same when evaluated in temporal gauge 
with a fixing function F Q [U] — J\ x 5[Uq(x), 1]. Writing 
the quantum mechanical trace over a complete set of 
states on the space H asTiO — ^ n {n\0\n) , one finds 

(0[U}) = (O[U]F [U]) 

= Z- lr ir lf Nt -* (fi(x)i(x,i)^ t (x)) ^ 

f* (fi(x)i(x,^t(x)) 1, (9) 

which has the spectral decomposition 

km {0[U}) = £ |<0| (n(x)i(x,i)nt( x )) |„0)|2 

JV t — >oo ' — ' \ /a/3 



xe 



-(E n -E )t 



(10) 



Here {|n )} is a complete set of eigenstates of T> and the 
matrix elements are assumed to be non-zero. 

On the other hand, one may also employ a fixing func- 
tion F L0 [U] = Yl x 5[n{x)Uo(x)n(x + 6),1], where now 
the composite Vb are brought to temporal gauge. This 
leads to the expression 

(0[U\) = (O[U]F L0 [U\) 



{f N *-*A(X, i)ap f| i(x, i)p a } , (11) 



where we have defined a modified transfer matrix Tlo 
with kernel 



K L o[U,U'] =e 



-Stiu,&n',u'] 



-S t [U u ' ,1,V n ] 



(12) 



and the notation CI' — Q [U']. According to the last equa- 
tion the integral kernels K ,Klo are simply related by 



K LO [U,U'}=Ko[U n \u' n "}. 



(13) 



One easily verifies that Klo shares most properties of 
Kq. In particular it is real, symmetric, K lo [U, U'] = 
Klq[U', U], and square integrable, so that Tlo is a self- 
adjoint, compact operator. It is gauge invariant under 
time independent tranformations g(x). Moreover, the 
positivity proof for To in [ pd| goes through unchanged 
for Tlq- Beyond the functional form of Kq, it only uses 
the fact that the product U'(x,i)U'(x,i) is in a funda- 
mental representation of SU(N), which is also true for 
C nt (x, i), U n ' (x,i). Thus Tlo has a strictly positive 
spectrum. Finally, the transfer matrix T acting on TL is 
again obtained by projection, T = PTlo- 

There then exists a complete set {|n i0 )} of eigenstates 
of Tlo to obtain the spectral decomposition of Eq. m~ 



km (O[c/])=^|(0|i Q/3 (x)|n i0 )|V 



(E'„-E )t 



(14) 



Eqs. (|T0|) ,(14) are two spectral representations of the 
same expectation value (0[f]) and of the form 



-E„t 



-E't 



(15) 



with a n . b n > 0, for every t. In the limit t — > oo only the 
n = 1 terms are retained. Expanding the exponentials 
one obtains in this limit 

^P«i(£ir = E { ^b 1 (E' 1 r. (i6) 



Both sides of this equation represent a power series in 
t of the same analytic function of t, hence the coeffi- 
cients of t m are the same for every m. It follows from 
the m = term that a\ — b\, while the first derivative 
of the equation with respect to t yields a\E\ = biE[, 
and hence E\ — E[. Going back to Eq. (^5|) and sub- 
tracting the n = 1 term on both sides, the procedure 



may be repeated to establish a n = b n and E n = E' n term 
by term. We thus conclude that Tlo has the same spec- 
trum as To , and hence the constructed correlator falls off 
exponentially with eigenvalues of the Hamiltonian Ho- 

The operator Eq. (|J) involves temporal Wilson lines 
which may be interpreted as propagators of static 
sources, and the field energies extracted from this correla- 
tor are those of a gluon field in the presence of sources. As 
the continuum limit is approached, the energies will thus 
diverge because of the well-known divergent self-energy 
contributions of the temporal Wilson lines, Fig. [l| In or- 
der to retain a finite continuum limit, the operator has to 
be modified such that no divergent mass renormalization 
is present. 



p-q 



FIG. 1. The Wilson line self-energy. Eq. ( fJ^ J enforces 
p — p — q = 0. 



This can be achieved by observing that the transfor- 
mation behaviour of Vq(x, t) in Eq. (||) is independent of 
the spatial coordinates. In the construction of the oper- 
ator Eq. (||) , instead of Vb we may then use its timeslice 
average 



Vb(t) = 5>o(x,t)/||^Vb(x,i)||, 



(17) 



which has been projected back into the group. The times- 
lice average corresponds to the Fourier transform of Vb 
with zero momentum. If this is done in every timeslice, 
the sources represented by the Vb cannot emit a gluon at 
one t and reabsorb it at some later t as in Fig.|l]. Hence, 
the mass renormalization of the static source is switched 
off, and the sources remain classical external fields. The 
presence of fields Vb then merely accounts for the trans- 
formation behaviour, but has no effect on the gauge field 
energies measured by the modified operator 



0[U] =Tr A l (x,0)Vb(0,i)A I (x,i)Vb t (0,t) 



(18) 



which has a spectral decomposition as in Eq. flU). The 
energies extracted from the expectation values of Eqs. (Q) 
and (|l8|) should then differ by a cut-off dependent shift 
due to the selfenergy contribution Fig. [I]. 

One may now ask to what extent these results depend 
on the particular choice of H[U]. Clearly, any 17 S SU(N) 



local in time and transforming as in Eq. (g) permits 
construction of the gauge invariant observable Eq. (18). 
From the spectral representation it follows that all such 
observables fall off with the same spectrum, while f2 only 
enters the matrix elements representing the overlap of 
the operator with the eigenstates. 

The construction of the composite link variable Eq. (j^) 
may also be viewed as fixing Laplacian gauge on each 
timeslice It is crucial that f2 depends only on spa- 
tial links to preserve the transfer matrix. In the language 
of gauge fixing this means that the gauge is incomplete, 
with a global factor h remaining unfixed between time- 
slices. It can be completed by imposing the further con- 
dition Vb(i) = 1, thus fixing h(t). In this particular gauge 
the operator Eq. (|lj) reduces to a gauge fixed gluon prop- 
agator, falling off exponentially with eigenvalues of the 
transfer matrix. Since the spectrum is unaffected by the 
particular construction of f2, the statement holds for all 
gauges employing a unique Q[U] local in time. For exam- 
ple, fixing the Coulomb gauge by the standard minimiza- 
tion of R[U] = J^xA 1 - 1/NTt (Up'*)] in every timeslice 
produces an with the desired properties and a residual 
freedom h(t) G SU(N). (Of course, this gauge condition 
has the problem that it does not determine Q uniquely 
||). On the other hand, Landau gauge is non-local in 
time and no positive transfer matrix is defined. 

Fourier transforming the spectral representation 



Eq. (14), one derives the Kallen-Lehmann representation 
in the limits N t — > oo, L d — > oo, 



(0[U\) 



P(E,P) = 



dE 

o 



P 



(2n) 



e- Et+ipx p(E,p), 



S(E~uj(p)) + p(E,p). 



(19) 



Hence the propagator in momentum space has a gauge in- 
variant pole defining a parton mass, m = cj(0) = Ei—Eq, 
with a residue \J Z(p) = (0|A(O)|p), while p(E,p) con- 
tains all higher states. 

It should be stressed here that this does not imply an 
asymptotically free colour charged state in the spectrum 
of the theory. All asymptotic one particle states satisfy 
Gauss' law without static charges, they are eigenstates 
of the projected transfer matrix T acting on the Hilbert 
space Ti. of gauge invariant functions. By contrast, the 
additional eigenstates \n L0 ) S Hq satisfy Gauss' law in 
the presence of static sources and are eigenstates of the 
transfer matrix To. 

All these considerations are easily extended to matter 
fields. In particular one may carry out the same analyis 
for the gauge invariant quark correlator 



Tr ( *(x, 0)fi(x, 0)Vb(0, i)fi f (x, i)*(x, t) 



(20) 

One would expect the mass extracted from its exponen- 
tial decay to coincide with the renormalization group in- 
variant quark mass computed in non-perturbative renor- 
malization schemes [p| . 



In summary, it has been shown that quark and gluon 
propagators exhibit an exponential decay governed by 
eigenvalues of the Hamiltonian in all gauges that are lo- 
cal in time, implying a corresponding gauge invariant 
and non-perturbative singularity structure in momentum 
space. In particular, the first excitation energy over the 
vacuum, viz. the pole closest to the origin, may serve as 
non-perturbative definition of a parton mass. This pro- 
motes propagator based definitions of Debye and mag- 
netic screening in the QCD plasma to a non-perturbative 
level, and opens similar questions for the zero tempera- 
ture case. The analysis may be extended to other Green 
functions and should reveal valuable information about 
the colour dynamics on any distance scale. Examples to- 
gether with detailed numerical experiments in three di- 
mensions will be reported elsewhere. 

I thank O.Bar, P. de Forcrand and U.-J. Wiese for nu- 
merous discussions, and P. van Baal for useful comments 
on the manuscript. 
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